Isothermal reentrant effect in a mesoscopic cylindrical structure of a superconductor 

coated with a normal metal layer 
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The coherent phenomena in mesoscopic cylindrical normal metal (N) - superconductor (S) struc- 
tures have been investigated theoretically. The magnetic moment (persistent current) of such a 
structure has been calculated numerically and (approximately) analytically. It is shown that the 
current in the N-layer corresponding to the free energy minimum is always diamagnetic. As the field 
increases, the magnetic moment (current) exhibits jumps at certain values of the trapped magnetic 
flux and the NS structure changes to a state with smaller absolute value of the diamagnetic moment. 
This occurs when the persistent current is unable to screen the external field. The magnetic moment 
increase stepwise and the system changes into a new stable state. The magnetic field penetrates 
into a larger volume of the N-layer. The new state has smaller absolute value of the diamagnetic 
moment. Experimentally, this is interpreted as the presence of a paramagnetic addition in the sys- 
tem (paramagnetic reentrant effect). The results obtained are in qualitative agreement with the 
experiments conducted by P. Visani, A. C.Mota, and A. PoUini, Phys. Rev. Lett. 65, 1514 (1990). 

PACS numbers: 74.45. -fC, 74.50. +r 



I. INTRODUCTION 

The Meissner effect in a mesoscopic cylindrical struc- 
ture consisting of a superconductor coated with a pure 
normal-metal layer has interesting features due to the 
coherent quantum effects in the normal metal. They are 
observable when there is a good contact between S and 
N constituents. This problem within the quasiclassical 
Eilenberger formalism was first investigated theoretically 
by Zaikin [11 . 

Recent advanced technologies of preparation of pure 
samples have enabled to investigate the coherence prop- 
erties of mesoscopic samples taking a proper account of 
the proximity effect [3] . The samples were superconduct- 
ing Nb wires with a radius R of tens of microns coated 
with a thin layer d of high-purity Cu, Ag or Au. Mota 
and co-workers [3| detected surprising behavior of the 
magnetic susceptibility x of a cylindrical NS structure 
(N and S are for the normal metal and the superconduc- 
tor, respectively) at very low temperatures ( T < 100 
mK ) in an external magnetic field parallel to the NS 
boundary. Most intriguingly, a decrease in the sample 
temperature below a certain point (at a fixed field) 
produced a paramagnetic reentrant effect: the decrease 
of magnetic susceptibility of the structure is changed to 
an unexpected grow. A similar behavior was observed in 
the isothermal reentrant effect in a field decreasing to a 
certain value below which the susceptibility started 
to grow sharply. It is emphasized in Ref. [4| that the de- 
tected magnetic response of the NS structure is similar 
to the properties of the persistent currents in mesoscopic 
normal rings. 
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There have been numerous attempts to explain the 
paramagnetic reentrant effect theoretically 0-13 • How- 
ever, the predicted amplitudes of the effect were too small 
( except for Q ) to account for the experimental facts. 
Fauchere, Belzig and Blatter [gj explain the large param- 
agnetic effect assuming strong repulsive electron-electron 
interaction in noble metals. The proximity effect in the 
N metal induces an order parameter will be shifted by 
TT from the order parameter A of the bulk superconduc- 
tor. This generates the paramagnetic instability of the 
Andreev states, and the density of states of the NS struc- 
ture exhibits a single peak near zero energy. The theory 
developed in Ref. |[6j essentially relies on the assumption 
of the repulsive electron interaction in the normal region. 

Maki and Haas 0] made the assumption that in no- 
ble metals (Ag, Au) p-wave superconductors may occur 
with a transition temperature of order 10 niK. Below Tc 
p -wave triplet superconductivity emerges around the pe- 
riphery of the cylinder. The diamagnetic current flowing 
in the periphery is compensated by a quantized param- 
agnetic current in the opposite direction thus providing a 
simple explanation for the reentrant effect. As in [^J, the 
authors of T\ also allow for paramagnetic current in the 
system, which flows in the opposite direction to the dia- 
magnetic current. Its amplitude is sufficient to explain 
the reentrant effect, but this theory says nothing about 
the temperature and field dependences of magnetic sus- 
ceptibility at ultra-low temperatures and in low magnetic 
fields. 

A theoretical basis for understanding the paramagnetic 
reentrant effect has been proposed in |8h1Q]. The theory 
is essentially based on the properties of the quantized 
levels of the NS structure. The Meissner effect is rather 
special in a superconducting cylinder coated with a pure 
normal- metal layer. The applied magnetic field gener- 
ates superconducting current in the surface layer whose 
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coincidence with the chemical potential of the metal. As 
a result, the state of the system suffer strong degeneracy, 
and the density of states on the energy of the NS sample 
experiences resonance spikes This contributes signif- 
icantly to the magnetic moment and causes a reentrant 
effect. Note that in j§] the calculation was performed for 
orbital susceptibility. In f6| the explanation involves the 
spin (Pauli) susceptibility of the system. 

In this study we calculated the free energy of the NS 
structure and its magnetic moment (current of magneti- 
zation) in the magnetic field. An approximate analytical 
calculation was supplemented with a numerical one based 
on the exact spectrum of Andreev levels in the NS 
contact. Our approach is not based on application of the 
Eilenberger equation. We calculate the thermodynamic 
potential to obtain the magnetic moment. 



II. THEORY AND RESULTS 



A. Spectrum of Quasiparticles of the NS Structure 



FIG. 1: Two classes of trajectories in the normal metal of 
the NS structure in a magnetic field: trajectories forming 
the Andreev levels (a): trajectories colliding only with the 
dielectric boundary (b ). This figure has been taken from j^. 



thickness is equal to the field penetration depth S. Si- 
multaneously, the Aharonov-Bohm effect generates per- 
sistent current (through the mechanism of the Andreev 
scattering of quasiparticles) in the normal layer near the 
NS boundary. If N and S metals are separated by a di- 
electric layer destroying the Andreev scattering, the addi- 
tional current disappears, and the Meissner effect returns 
to its usual form. The levels with energies no more than 
A (2A is the gap of the superconductor) appear inside 
the normal metal bounded by the dielectric (vacuum) on 
one side and the superconductor on the other side. Be- 
cause of the Aharonov-Bohm effect [ll| , the spectrum of 
the NS structure in a weak field is a function of the mag- 
netic flux. The specific feature of the Andreev quantum 
levels of the structure is that by varying field H (or tem- 
perature T) each level in the well periodically comes into 



Consider a superconducting cylinder with the radius 
R which is concentrically embedded with a thin layer d 
of a pure normal metal. The structure is placed in a 
magnetic field (0,0,11) oriented along the symmetry 
axis of the structure. It is assumed that the field is weak 
to the extent that the effect of twisting of quasiparticle 
trajectories becomes negligible. It actually reduces to the 
Aharonov-Bohm effect i.e., allows for the increment 
in the phase of the wave function of the quasiparticle 
moving along its trajectory in the vector potential field. 

We proceed with a simplified model of NS structure in 
which the order parameter magnitude changes stepwise 
at the NS boundary ( A (x) = 6(-a;), Im (A) = 0, a 
bulk superconductor in the region a: < and a normal 
metal layer in the region < x < d ). It is also assumed 
that the magnetic field does not penetrate into the su- 
perconductor. The coherent properties observed in the 
pure normal metal can be attributed to its large "coher- 
ence" length (r) = :^^ffff ( is the Fermi velocity, 
ks is the Boltzmann constant) at very low temperatures 
(Cat {T) >• d). Besides, the spectrum of quasiparticles 
was obtained assuming a negligible curvature of the NS 
boundary. 

One can easily distinguish two classes of trajectories, 
inside the normal metal. One of them includes the trajec- 
tories which collide in succession with the dielectric and 
NS boundaries (Fig.l). The quasiparticles moving along 
these trajectories have energies \E\ < A and are local- 
ized inside the potential well bounded by a high dielectric 
barrier ( ~ leV ) on one side and by the superconducting 
gap A on the other side (A = 3.56- kBTj2, A(Nb) « 
1.42 meV). On its collisions, the quasiparticle is reflected 
specularly from the dielectric and experiences the An- 
dreev scattering at the NS boundary [ij]. We introduce 
an angle a at which the quasiparticle hits the dielectric 
boundary. The angle is measured in the positive direc- 
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FIG. 2: The dependence of the density of states of the NS 
structure on the magnetic flux $ ( i? = Ep = ). Normal- 
ization was performed for the flux ^max corresponding to the 
highest value of v ($) ( $,nax ~ 2.175 ). 



tion from the normal to the boundary (Fig.l). In this 
case the first class contains the trajectories with a vary- 
ing within the range — Qc < a < ac { ac is the angle at 
which the trajectory touches the NS boundary, sin(ac) 



= ■j^:^ )• Another class includes the trajectories whose 
spectra are formed by collisions with the dielectric only, 
i.e., the trajectories with a > ac- The two groups of tra- 
jectories produce significantly different spectra of quasi- 
particles. The distinctions are particularly obvious in 
the presence of the magnetic field. The trajectories with 
a < Qfc form a spectrum of Andreev levels which con- 
tains an integral of the vector potential field. The spec- 
trum characterizes the magnetic flux through the area of 
the triangle between the quasiparticle trajectory and the 
part of the NS boundary. It also determines the mag- 
nitude of the screening current produced by "particles" 
and "holes" in the N layer. These states are responsible 
for the reentrant effect. The trajectories with a > Qc do 
not collide with the NS boundary. The states induced by 
these trajectories are practically similar to the "whisper- 
ing gallery" type of states appearing in the cross section 
of a solid normal cylinder in a weak magnetic field Til], 
The size of the caustic of these trajectories is of the 
order of the cylinder radius, i.e., they correspond in high 
magnetic quantum numbers. The spectrum thus formed 
carries no information about the parameters of the su- 
perconductor, and it is impossible to meet the resonance 
condition in this case. These states make a paramagnetic 
contribution in the thermodynamics of the NS structure 
but their amplitude is small ( ^ 1/ {kp ■ R) )■ It is there- 
fore discarded from further consideration. Our interest 
will be concentrated on the trajectories with \a\ < ac- 

The spectrum of quasiparticles of the NS structure can 
be obtained easily using the multidimensional quasiclas- 
sical method generalized for the case of the Andreev scat- 
tering in the system [15| , [l6( . After collision with the NS 
boundary the "particle" transforms into a "hole". The 
"hole" travels practically along the path of the "particle" 
but in the reverse direction (Fig.l). 

The spectrum was derived by quantizing the adia- 
batic invariant § P ■ ds , where P = p+ ^A, A = 
{0,Ay {x),0), Po = po- for a "particle" and Pi = 

Pi+ —A for a "hole". Note that each collision with the 
NS boundary multiplies the wave function amplitude of 
the quasiparticle by a factor of exp (— iarccos (i?/A)). 
Let £o be the length of the quasiparticle trajectory be- 
tween the collisions at the boundaries of the N layer. We 
thus [1], [§| arrive at the expression for the spectrum of 
the Andreev levels in the NS structure: 



En (g, a, $) 



Co 



En (g, a, <&) \ tan (a) 



(1) 



Here Vi(q) = -^/pir^^^/m*, Cq is length of the quasi- 
particle trajectory, pf is the Fermi momentum, q is the 
quasiparticle momentum component along the cylinder 



axis ( \q\ < pf ), rri* is the effective mass of the quasipar- 
ticle, and $0 = hc/2e is the superconducting flux quan- 
tum. The factor $ appearing in the last term in Eq. ([1]) 
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Ay (x) dx 



(2) 



which is dependent on the vector potential field A = 
{0,Ay{x),0). The spectrum of Eq. ([T]) is similar to 
Kulik's spectrum [l3| for the current state of the SNS 
contact. However, Eq. ([T]) includes an angle-dependent 
magnetic flux instead of the phase difference of the con- 
tacting superconductors. 

The length of the quasiparticle trajectory ( 2AB ) is 
readily found from Fig.l using the sine and cosine theo- 
rems: 



(cos (a) — \/sin (ar)'^ — sin (a)^ \ 
— -4^. — (3) 
1 - sin (ac) j ^ ' 

where sin(ac) = < a < Uc ■ The spec- 

trum in Eq. ([T]) was derived assuming that the mean 
free path of the quasiparticles was much longer than the 
cross-section perimeter of the cylinder and the require- 
ment d <^ R was obeyed. In this limit £o — 2AB ~ 



2c?/ cos (a) lim AB = d/ cos (a)] , i.e., the radius 

R drops out from the expression for the spectrum. Al- 
though the boundary curvature of the sample is disre- 
garded, the information about its cylindrical geometry 
is retained through a correct choice of the limits of in- 
tegration for the angle a: —ac < a < etc- Putting Cq 

FIG. 3: Free energy (normalized per value Q. (4>max)) as a — 2d/ cos (a), we obtain the following expression for the 

function of the flux $. spectrum ( as in [9| ): 



(g,a,$) 



7rfi,VL((?) cos (a) 
2d 



1 ( E^[q,a,^) 

n H — arccos 

TT V A 



tan (a) 



(4) 



r 



The spectrum in Eq. Q has an important feature 
[1], [i]. As the "phase" $ Eq. ^ changes, the density 
of states exhibits resonance spikes. Every time when the 
Andreev level coincides with the chemical potential of the 
metal, the state of the NS structure suffers strong degen- 
eracy showing up spikes. The dependence of the density 
of states upon the magnetic flux calculated numerically 
for the NS system is illustrated in Fig. 2. 

Note that in [3] the diamagnetic current of NS 
structure was calculated using Uc = 7r/2 instead of the 



upper limit of integration for the angle ( ac < 7r/2 ) and 
assuming implicitly an infinitely large number of Andreev 
levels. Therefore, these results cannot be conformed with 
the experimental findings. The reason is not only that 
the calculation was made for a flat geometry rather than 
for a curved NS boundary. Numerical analysis shows that 
an adequate interpretation of the experimental magnetic 
moment-field dependence is possible only with a proper 
choice of the upper limit of integration with respect to 
a. If a > Qfc or ckc = 7r/2, the consideration includes 
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effectively the states unrelated to the Andreev levels. 



B. Self-consistent equation 

To calculate the "phase" $(r,H) from Eq. we 
should know the distribution of the vector potential field 
inside the normal metal. Zaikin [l| has shown that the 
proximity effect caused the Meissner effect leads to an 
inhomogeneous distribution of the vector potential field 
over the N layer of the structure: 



Ay (x) — ^qRx + fj,Qj ■ x(d— x/2) 



(5) 



where /ig is the permeability of free space ( the SI sys- 
tem of units is employed, the geometry of the proxim- 
ity model system is the same as in [19] Fig 1). This 
expression can be obtained from the Maxwell equa- 
tion rot [rot (^A^^ — ji^j = assuming that 
the current density is uniform over the cross-section of 



the conductor and the boundary condition A 



rot 



x=0 



0, 



= i3(0,0,/iQH) is met. The fact that the 

current density is constant in the N-layer follows from 
spatial homogeneity of the density of Andreev levels over 
the whole thickness of the N-layer. In cylindrical geome- 
tries, if the N-layer thickness is not thin compared to the 
radius {ji ^ R ), the current density is not constant in 
space pQl]. 

The magnetic moment per unit length of the N-layer 
M (H) ~ dn(TMiiy) ^ 2_component ) and the current 

density j are related via a relation 



M (T, H) 



•x ]{r) 



dV ^ 



Vn 



-I dn 



(6) 



where Vn is a volume of the N-layer unity height, il (T, <i>) 
is the free energy per unit length. Then, the current 
appears from Eq. ([S]) is a function of the magnetic flux 
$ and temperature: 



1 



dn (r,$(H)) 



nR'^d- 



Mo 



dR 



(7) 



We can write down the self-consistent equation for 
$(T,H) using Eqs. ©, ® and Q: 



$ (T, h) = h + 7] • M* ($) 



9$ (T, h) 



where h 

dn 



H 



_ <S>o 



— ^ [Ml- To describe the field effect on the magnetic 
moment 



ah 

rf' M* ($) 



(8) 



M (T, H) 



M* (T, a$ (T, h) 



dh 



(9) 



of a NS structure, it is necessary to find the dependence 
$ (T, h) from Eq. ^ . After calculating the free energy 
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FIG. 4: The magnitude of M* ($) / |M* ($„ 
of the flux $. 



as a function 



57 (T, <i>) from the spectrum of Eq. Q , we can estimate 
the magnetic moment Eq. ^ using a solution of the dif- 
ferential equation ©: M (T, H) = (T, h) - h). 



In this article, we used the "thermodynamic" approach 
Eq. ([7]) which leads to the first order differential equation 
([5]) for the function <I> (T, h). However, another approach 
based on of the Eilenberger formalism ( [l^ ) yields 
an algebraic self-consistent equation ( Eq. 22 in (fl^) 
) for the "phase" $(r,h): 



$ (T, h) = h + const ■ j ($) 



(10) 



( in notation Eq., ([5]) ). In this equation the function 
j ($) is described by the expression Eq., 13 in [l|. Clearly, 
both approaches (HI), (fTO)) will lead to quite different de- 
pendencies of M (T, H) on magnetic field and tempera- 
ture. In our point of view, the self-consistent equation 
(fTU)) cannot be applied to the cylindrical NS structures. 
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While deriving the expression for a current the au- 
thor assumed that ac = 7r/2 (the case of a plane). This 
means account for the contribution non-Andreev states 
( a > ttc). For the calculation of thermodynamic po- 
tential in equation (jS| we shall use the actual magnitude 
of the parameter ac ( sin{ac) = R/{R + d) ). The ap- 
proximation d « R was used by us in the derivation 
of the spectrum Q only. In other words, the neglect of 
the curvature of cylindrical samples (i.e. the path length 
of quasiparticles was choosen as d/cos(a) (jS])) does not 
entail the need of account for the states with a > for 
the cylindrical NS of structures (Fig.l). 



where the summation is over the spin variable s = ±1 
and all the states related to the quasiparticles trajectories 
with |q;| < cic in Eq. Then, we obtain the following 
expression for the free energy per unit length L 



C. Analytical estimation of the magnetic moment 
of the NS structure 

We proceed from the expression for the free energy of 
a NS contact: 

ln(l + e "B-T \ (11) 



R-keT 



n— — oo 



PF 



In ( 1 + exp ( — 



^„ (g,a,^) 



Pp — cos (a) dqda 



(12) 



where the energy En (q, a, $) is given by the exact ex- 
pression for the spectrum in Eq. (jj]). For simplification, 

we introduce the dimensionless quantities = , 

— 1 < £« < 1, ( A is the superconduct- 



2d-A-m* ' 



ing gap ) and perform the change of variables {q, a} 
{u,v}: 



1- ( ^ 

PF 



■ cos (a) 



V — (J 



'l-(^) .sin (a) 

The spectrum and the free energy become 

e„ = [nil + arccos (£«)] • u — $ • w. 



r!($) = cir5] 



dn ^2 cosh ^2cjr)) dudv 



(13) 



(14) 



(15) 



where ci ^ _2^^ • {^)\ C2 = ^, < < a, 

— (7 sin (ofc) <v< cr sin (ac), — £„ (w, <&), an inte- 
gration domain S" is a sector of a circle of radius a . In the 
expression (|15|) we also took into account the symmetry 
of the spectrum in Eq. ([1] 



(16) 



-|n| (w, V, $) = -£|„|_i (u, -U, $) . 



Making use of the relation 
d£n 



(17) 



we evaluate the derivative of the free energy with respect 
to the flux M* ($) = 



M* ($) - C3 ^ 



ri=0 



'tanh (2^) v^l - eidudv 



(18) 

where C3 = ^ (f^)^ Eqs. ®, © and dH]) fully de- 
termine the non-linear magnetic response of a cylindrical 
NS structure to an externally applied magnetic field H. 

The integral expression of Eq. suggests that 

M* ($) is the odd function of the flux $: M* ($) = 
—M* (— $). A hnear term of the function M* ($) has 
been determined from an approximate estimation of the 
integral in Eq. (fT5|) . This calculation is similar to that 
in the Attachment of [1]. The final expression for the 
magnetic moment is 

, , ^ In cosh (^) , , 

M (T, h) ~ Mo — ^ (19) 

TT-n 
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where n 

Mo - - 



C3 • CT" 



Ta 

(ft) 



2TTd-k 



is Andreev temperature, 
$(r,h) • ^5^, the "phase" 

$ (T, h) is a solution of the differential equation (|8]), uq 
is the number of Andreev levels in the potential well ( 

no ~ 5^-^^ tan (ckc) )• Eq- (fT9)) shows that the mag- 
netic moment is diamagnetic in the range of small fields 
( $ (h) = const ■ h, const > ) and allows for the contri- 
butions of " particles" and " holes" . 



D. Numerical results 

Let us compare two approaches described above for the 
calculation of the magnetic moment of the NS structure. 
Fig. 5 shows the function M* ($) and dependency of the 
current on the magnetic flux obtained in the Green's 
function approach. For comparison, we obtained the de- 
pendence AI* ($) at the same value ac = 7r/2 as was 
used in the derivation of the formula j ($) in [3]. In 
the initial part (linear in the $) both curves coincide. In 
this approximation ($ (h) = const ■ h) the self-consistent 
equation © turns into Eq. PU)) . Thus, at small values 
of the magnetic field we would obtain the same field de- 
pendence of the magnetic moment M (T, h) for the NS 
structure in both approaches. However, in large fields 
the behavior M {T,h) is quite different. To calculate 
M (T, h) from Eq. © , we have used the following physi- 
cal values of the NS structure: R = 8.3 ^m, d — 3.2 ^m, 
{ac = 36°), v^. (Au) = 1.4-10^ cm/s, A (Nb) = 1.42 meV 
( a = 0.644, r/ • C3 = 5.3 • 10^, Hq = 51A/m = 0.64 Oe 
). The selectedparameters are close to those used in the 
experiment [1, 13, [H]. 

The results of calculation according to formulas ([15]) 
and (dH]) are illustrated in Fig. 3, 4. While plotting 
Fig. 3, the nonzero quantity 17 ($ = 0) was omitted. 
The dependence M* ($) = — 35 (Fig. 4) crosses the ab- 
scissa thereby determining singular points of the differ- 
ential equation ([5]) . The dependence <i> (h, T) calculated 
through numerical solution of the self-consistency equa- 
tion ^ exhibits jumps and is illustrated in Fig. 6a for 
the branches corresponding to the minimum of the Gibbs 
free energy (23| : 



1 

2/io 



Vn 



dV (20) 



where B — rot (1), H = i?(0,0,H). The magnetic mo- 
ment M (h) and the free energy ft (h) as functions of the 
magnetic field are shown in Fig. 7a, Fig. 6b. Each 
jump Ail of the free energy (see Fig. 6b) is accompa- 
nied by the jump of the magnetic moment AM (see Fig. 
7a) in such a way that the Gibbs free energy (f^Ol) is a 
continuous function of versus the magnetic field h. We 
have not performed an analysis of behavior Gibbs free 
energy (|20p near the points where the magnetic moment 
has jumps because it is beyond the semi-classical approx- 
imation adopted in this article. 
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FIG. 5: The magnitudes Af* (<1>) / |M* (0.1)| and 
J (<I>) / |jf (0.1)1 as the functions of the flux $ (for expla- 
nation see text). 



III. CONCLUSIONS 



The goal of our study was to interpret the experiments 
performed by A. C. Mota et al., who detected an 

anomalous behavior of the magnetic susceptibility of the 
NS structure in a weak magnetic field at millikelvin tem- 
peratures. Previously [^,9], the anomalous behavior of 
the NS structure was attributed to the properties of the 
quantized Andreev levels depending on the magnetic flux 
that varies with the temperature and magnetic field. We 
used the "thermodynamic" approach for the calculation 
of the magnetic moment of normal region of NS structure. 
Within the framework of the self-consistent equation (|8| , 
we have managed to trace the role of the parameter Uc 
in thermodynamics of NS of structures (Fig. 4 and Fig. 
5). Failure connected with the use of the quasi-classical 
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FIG. 6: a) the dependence of <1> (T, h) on the magnetic field h. 
b) the relative free energy (h) /SI (hmax) (h max J — ^max J 
as a function of the magnetic field h. 



FIG. 7: a) The magnetic moment of the NS structure versus 
the magnetic field h. b) Isothermal dc-magnetization curve 
at T = 7 mK for the sample 41AuNb Qf^. Used by courtesy 
of A. C. Mota. 



Green-function technique [H] for the explanation of the 
experimental data (Mota et al) is a consequence of ac- 
count for states which are non-Andreev ones ( a > etc 
) for the cylindrical NS structures. Geometrically, this 
can be seen from figure Fig. 1. The quasiparticle tra- 
jectories ( a > ttc ) hitting the dielectric boundary only 
are responsible for the paramagnetic current of a small 
amplitude ( ^ l/(fcF • -R)) (we neglected this current). 

The proximity effect is crucial for the reentrant effect. 
The amplitude of the resonance spikes in the density of 
states strongly depends on the probability of the Andreev 
reflection at the NS boundary. It is therefore assumed 
that the normal metal and the superconductor are in a 
good electric contact. The spectrum Eq. (j4|) was ob- 
tained by the method of multidimensional quasi-classical 
approach [l^, [l^. In doing so, we assumed (i) the condi- 



tion of smallness of N-layer thickness in comparison with 
the radius of the cylindrical superconductor, (ii) the va- 
lidity of the model of a stepwise- varing order parameter of 
the structure, (iii) the independence of A on the magnetic 
field. This permitted us to pass over from the curved NS 
boundary to a flat one. The information about the cylin- 
drical geometry of the sample was retained because the 
critical angle at which a quasiparticle hits the dielectric 
boundary a is smaller than olc- The problem was further 
simplified by assuming that the reflected quasiparticle 
performed a reciprocating motion, i.e., a "particle" and 
a "hole" pass along the same trajectory but in opposite 
direction. Actually, here there exists a lot of quasi recip- 
rocating trajectories with the energies near e ^ 0. These 
trajectories lead to the spikes in the density of states and 
were taken into account for numerical computation. The 
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numerical calculation shows that the nonlinearity of flux 
<i>(T,h)- field dependence ( T= constant ) (Fig.6a) 
gives rise to quite interesting features of M (T, H) . The 
magnetic moment in the N layer appears to be always dia- 
magnetic. The paramagnetic contribution to current (the 
paramagnetic reentrant effect) was not detected. How- 
ever, we have obtained the stepwise change in absolute 
value of the magnetic moment with increasing magnetic 
field (Fig. 7a). This behavior can be interpreted as an 
appearance in the magnetic moment the paramagnetic 
additives. A behaviour of the NS structure changes from 
one stable state to another and the magnetic field pene- 
trates further into a bulk of the N layer. The new state 
has smaller absolute value of the diamagnetic moment, 
which is interpreted experimentally as an evidence of a 
paramagnetic addition in the system (Fig. 7b). When 
the field grows further, |A/(H)| increases again until its 
new value makes the system to jump to a next stable 
state with a smaller absolute value of the diamagnetic 
moment, and the magnetic field penetrates deeper inside 
the normal metal (Fig. 7a). The number of the moment 
jumps depends on the number of Andreev levels in the 
NS structure. Under the isothermal condition, the values 
of the magnetic field at which jumps occur, do not coin- 
cide when the magnetic field changes from small to larger 



values and in the opposite direction because of a special 
dependence of the Gibbs free energy on the field. This 
sort of hysteresis was observed experimentally in 0, [23 | 
(Fig.7b). 

Numerical comparison between data presented at 
Fig.7b and Fig. 7a shows that M(r,H) (Fig. 7a) 
gives the qualitative description of the experimental data 
which obey the scaling rule: H^^'P/Hf ^ ~ ^caic jY^caic ^ 
5/2. For the quantitative description of the temperature 
and magnetic field dependence of the magnetic moment 
NS of structures, it will be important to take into account 
the exact spectrum of Andreev levels, the latter is sup- 
posedly possible within the framework of the Bogoliubov 
- de Gennes equations only. 

Note that our consideration was entirely based on 
the model of free electrons without account of strong 
electron-electron repulsion. 
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